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NOMENCLATURE 


Cross-sectional  area  of  the  bar 

Phase  velocity,  c  =  sec  5/2,  in./sec 


Constant 

Integration  constants 


Complex  modulus,  psi 

Real  part  of  complex  modulus,  psi 


Imaginary  part  of  complex  modulus,  psi 
Magnitude  of  complex  modulus,  psi 
Frequency  of  vibration  when  Re  =  0,  He 

<pr 

Imaginary  part  of  the  ratio  of  the  acceleration  of  the  driven  end  of  the  bar  to  that 
of  the  (tee  end 

Length  of  the  test  bar,  in. 

End  mass,  lb-sec 2/in. 

Mode  of  vibration 

-  ^1  -  i  tan  ^  j ;  see  also  definition  below  eq  4 

Absolute  value  of  the  ratio  of  the  acceleration  of  the  free  end  of  the  bar  to  that  of 
the  driven  end 

Measured  acceleration  ratio  when  Re  =  0 
Mass  ratio 

Real  part  of  the  ratio  of  the  acceleration  of  the  driven  end  of  the  bar  to  that  of  the 
free  end 

Time,  sec 

Displacement  at  any  section  of  the  bar  as  measured  on  the  x  -  y  coordinate 
system,  in. 

Amplitude  of  displacement  u,  in. 

Displacement  at  the  fixed  end  of  the  bar,  in. 


Amplitude  of  displacement  u0,  in. 

Axial  coordinate,  in. 

Normal  coordinate,  in. 


NOMENCLATURE  (Cont’d) 


moj2 

pAE(ii j) 


) 


Angle  by  which  strain  lags  stress,  radians 
Strain,  in. /in. 

Amplitude  of  strain  in  a  sinusoidal  excitation,  in. /in. 

Frequency  ratio,  £  =  <oL/c 

Frequency  ratio  when  Re  =  0 

Mass  density,  lb-sec 2/in. 4 

Axial  stress,  psi 

Amplitude  of  stress  in  a  sinusoidal  excitation,  psi 

Phase  angle  between  bar  end  absolute  displacements,  radians 

Exciting  angular  frequency,  rad/sec 

Exciting  angular  frequency  when  Re  =  0,  rad/sec 
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INTRODUCTION 

The  U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory  (USA  CRREL)  employs  a 
test  technique  for  determining  the  complex  moduli  and  damping  of  frozen  and  non  frozen  soils  under 
vibratory  loads.  This  involves  submitting  an  upright  cylinder  of  the  material  to  vibration  at  the 
lower  end  with  the  upper  end  free.  Input  and  output  wave  characteristics  are  measired  by 
accelerometers  fastened  to  a  base  plate  and  top  plate,  respectively.  Other  investigators  are  known 
to  employ  similar  techniques  for  testing  a  variety  of  materials.  In  the  analysis  of  test  measurements 
to  obtain  the  desired  properties  of  the  material  the  authors  have  determined  the  effect  of  the  top 
end  plate.  They  show  that  the  mass  of  the  end  plate  in  comparison  with  the  mass  of  the  sample  has 
a  significant  effect  on  the  measured  moduli  and  damping  properties  of  the  material. 

A  convenient  method  of  measuring  the  complex  modulus  of  a  linear  visooelastic  material  over 
the  audiofrequency  spectrum  is  to  apply  a  harmonic  displacement  to  one  end  of  a  bar  of  the  material 
and  measure  the  ratio  of  end  accelerations.  The  problem  has  been  considered  by  Lee  (1963)  and 
Brown  and  Selway  (1964)  whose  orientation  was  directed  to  materials  as  diverse  as  soils  and  poly¬ 
mers.  The  solutions  given  by  these  authors  specify  a  free-end  boundary  condition.  However,  many 
experimenters  using  this  technique  have  found  it  convenient  to  measure  the  end  displacement  with 
an  accelerometer.  The  work  presented  here  accounts  for  this  end-mass  effect  and  indicates  the 
deviations  one  may  expect  from  the  simpler  free-end  theory. 

The  theoretical  work  presented  here  is  supplemented  by  experimental  results  which  indicate  the 
applicability  of  the  theory. 


THEORY 

Derivation  of  eqaatioss  for  displacement,  strata  aad  stress 

The  equation  describing  the  motion  is  most  easily  obtained  by  assuming  the  x  and  y  axes 'fixed 
in  the  bar  (see  Fig.  1)  at  the  driven  end  x  =  0  and  the  system  given  a  displacement  u0  *  U0 
exp  (ieot). 


The  equation  of  motion  is 
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where  o  is  the  uniaxial  stress,  p  is  the  mass  density  and  u  is  the  axial  displacement  of  a  point  in 
the  bar  measured  relative  to  the  moving  coordinate  system.  Taking  the  stress  at  a  point  in  the  bar 
as  a  =  &  .exp  (ifc*t)  and  the  strain  as  t  =  T  exp  [i(wt  -  5)),  the  constitutive  law  may  be  written  as 


~  =  Z  exp  (iS)  =  E *  exp  (iS)  =  E(ia>).  (2) 
Taking  u  =  ff  exp  (i&»t)  and  eq  2  in  the  form 
a  =  (3) 

eq  1  goes  over  into  the  ordinary  differential  equation 

~  +  p2*  =  -p2u0 

dxK 

where  ps  «  pw2/£(/<u). 

The  solution  to  eq  4  is 


Figure  1.  Coordinate  system. 


(4) 


+  uo  =  Cjoospx  +  CgSinpx 

where  C  j  and  C2  are  obtained  from  the  boundary  conditions 

u(0,t)  =  0 


d* 

Ao(L,t)  =  -m — (u  +  un)__,. 

dt*  L 


(6) 


A  is  the  cross-sectional  area  of  the  bar  and  m  is  the  end  mass.  Applying  these  boundary  conditions 
to  eq  5  the  displacement  solution  is 


tKx.to) 


008  PX  + 


(tanpL  +  y  ^ 
1  -  y  tanpL  / 


sin  px  -  1 


(7) 


where 


m«8 

Y  ~  pAE(iw)  ' 

The  stress  and  strain  at  any  point  in  the  bar  are  respectively 


(8) 


a  =  E*  exp (iS)( 


(9) 
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and 


(  =  j  cos  px  -  Bln  pxj  exp  (lot)..  (10) 

The  solutions  in  eq  7*10  may  be  put  in  more  meaningful  form  by  substituting  the  expressions  for  p 
and  y  and  separating  tbe  right-hand  side  into  real  and  imaginary  parts.  The  complex  result  repre¬ 
sents  the  magnitude  of  the  displacement,  strain  or  stress  and  the  phase  relative  to  the  base  dis¬ 
placement  at  any  point  x. 

Solution  is  terns  of  the  bar  aid  aceelwstloa  ratio 

A  simple  relationship  may  be  found  for  the  ratio  of  bar  end  displacements  (or  accelerations) 
that  is  useful  in  experimental  measurement  of  the  oomplex  modulus.  Rewriting  eq  7  for  x  »  L  gives 


n (L,<o)  +  u0 

sec  pL 

1  -  y  tan  pL 

(ID 


It  is  convenient  to  define  the  frequency  ratio 


£  = 


<dL 

c 


(12) 


where  c  is  the  phase  velocity  y/E*/p  sec  5/2.  Using  eq  8  and  12  and  the  definition  of  p,  eq  11  may 
be  put  in  the  form  of  a  real  and  imaginary  part 


Ur 


U{L,u»)  +  UQ 


Re  +  ilm 


where  it  may  be  shown  after  some  algebra  that 
Re  =  cosh 

and 


((Mi) 


(cos  £  -  R  £  sin  £)  +  R  £  tan  -  cos  £  sinh 


Im  =  sinh 


M) 


(sin  £  +  R  £  oos  £)  +  R  £  tan  -  sin  £  cosh 


(f“i) 


(13) 


(14) 


(15) 


where  R  is  the  mass  ratio,  m/pAL.  Equations  11  and  15  are  also  valid  for  large  values  of  8,  no 
simplifying  assumptions  having  been  made. 


Use  of  end  acceleration  ratio  to  sMasure  the  complex  modulus 

Equations  14  and  15  suggest  an  experimental  technique  to  measure  the  complex  modulus  of  a 
linear  viscoelastic  material.  The  experimental  technique  for  measurement  of  in-phase  and  quadra¬ 
ture  components  of  the  response  is  within  the  state  of  the  art  with  commercially  available  equipment. 
Measurement  of  the  complex  end  displacement  ratio  (or  equivalently,  the  acceleration  ratio)  yields 
experimental  values  for  Re  and  Im.  Substitution  of  these  two  values  in  eq  14  and  15  yields  two 
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simultaneous  transcendental  equations  which  may  be  solved  numerically  for  the  two  unknowns  £ 
and  tan  S/2  for  any  mass  ratio  R.  Having  then  solved  for  and  tan  S/2,  the  complex  modulus  may 
be  easily  obtained  from  eq  12  and  the  definition  of  the  phase  velocity  c.  Specifically 

E*  =  pc2  cos2  |  =  g-  cos  g j  *  (1®) 

Hence 

=  £*  cos  5  (17) 

and 

£2  =  £*  sin  5  (18) 

where  £j  and  £ 2  are  real  and  imaginary  parts  of  the  complex  modulus. 

Measurement  of  the  complex  modulus  at  90°  phase  shift 

A  simple  experimental  method  to  determine  the  complex  modulus  in  the  vicinity  of  the  bar 
resonances  was  suggested  by  Lee  (1963)  and  Brown  and  Selway  (1964).  The  phase  relationship 
between  the  end  displacements  is  given  by  eq  13.  If  4>  is  the  angle  between  the  displacement  of 
the  driven  end  to  the  free  end  of  the  bar 


<t> 


(19) 


It  follows  that  when  Re  =  0  there  is  a  90°  phase  shift  which  is  easily  measured  experimentally 
without  sophisticated  equipment.  For  this  case  eq  14  and  15  reduce  to 


Re  =  cosh 


(f“l) 


(cos  -  R£'  sin  £')  +  Rf'  tan  -  cos  <f'  sinh 

6 


(f,an!) 


and 


Im  =  sinh 


(fan!) 


(sin<f'  +•  R<f'cos£')  +  R  <f'  tan  -  sin  cosh 

2 


1 

\  2/ 

Q' 

0  (20) 


.  (21) 


Q '  is  the  measured  acceleration  ratio  when  Re  0;  the  frequency  ratio  at  this  point  is  defined  as 
f  and  the  frequency  as 

The  experimental  procedure  is  to  adjust  the  frequency  until  the  phase  relationship  is  90°;  at 
this  point  the  frequency  and  the  acceleration  ratio  Q'  are  measured.  Using  Q',  eq  20  and  21  may 
be  solved  numerically  for  and  tan  5/2;  hence  the  complex  modulus  may  be  calculated  using 
eq  16-18.  This  method  limits  the  data  to  a  specific  frequency  in  the  vicinity  of  the  bar’s  resonant 
frequency.  does  not  generally  coincide  with  £  at  resonance  as  is  shown  later  in  this  report. 

A  computer  program  or  a  set  of  curves,  both  given  in  this  report,  may  be  used  to  solve  for 
and  tan  5/2  using  experimental  data.  A  computer  program  in  Fortran  IV  is  given  in  Appendix  A  to 
solve  eq  20  and  21.  This  program  uses  the  Newton-Raphson  method  (see  Scarborough,  1955)  to 
solve  for  and  tan  5/2.  The  program  reads  R,  Q\  p,  L,  mode  and  base  amplitude  and  prints  out 


LONGITUDINAL  FORCED  VIBRATION  OF  VISCOELASTIC  BARS  WITH  END  MASS 


5 


c,  £',  tan  5/2,  E*,  E E 2  and  base  stress.  Alternatively,  curves  are  given  in  the  section.  Com¬ 
puter  generated  curves  and  discussion  (p.  7),  to  obtain  <f‘  and  tan  5/2  directly  using  experimental 
values  of  R  and  Q'. 

Stress,  strain  and  displacement  for  any  value  of  x 

The  separation  of  eq  7-10  into  real  and  imaginary  parts  to  obtain  useful  formulas  for  stress, 
strain  and  displacement  as  a  function  of  x  leads  to  complicated  algebraic  expressions.  This  effort 
may  be  circumvented  by  making  use  of  the  computer’s  ability  to  do  complex  arithmetic  directly. 

The  stress,  strain  and  displacement  were  evaluated  at  five  stations  along  the  bar  for  various  mass 
ratios.  Typical  results  are  given  in  Table  I.  The  computer  program  is  given  in  Appendix  B.  Values 
given  are  for  a  soil  with  E*  =  78,000  psi  and  tan  5/2  =  0.06. 

The  maximum  stress  occurs  at  x  =  0  for  the  first  three  modes.  An  expression  for  this  stress 
derived  from  eq  9  and  10  is 


<nO,w)  =  E(icj)U0p(  t*npL  + y). 

\  1  -  ytan  pL  / 


(22) 


Applying  the  complex  definitions  of  p,  £(/«)  and  y,  and  considering  the  simple  case  when  the  90° 
phase  shift  occurs  (Re  =  0),  one  has 


<r{0,a>)  =  V”[run  |  -  IT  -  ftan|  +  /J J 


there 

R  =  (cob  -  R£‘  sin  £')  sinh  tan  ^  j 

T  =  (sin  +  R£'  cos  £')  cosh 


+  Rf'  tan  -  cos  cosh 


(f,“n0 


+  /?£'  tan  -  sin  sinh 
*5 


(f“l) 


(23) 


and  Q'  is  the  experimentally  measured  value  of  the  acceleration  ratio  at  a  90°  phase  shift.  The 
magnitude  of  the  stress  is  given  by 


<7(0, w)]  =  Re2(<7)  +  Im2(<r)  . 


(24) 


This  calculation  is  built  into  the  standard  program  for  calculating  tan  5/2  and  E*  from  experimental 
data  (see  Appendix  A). 
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Redaction  of  theory  to  earlier  work 

For  tan  5/2  =  0  eq  14  and  15  reduce  to  the  equation  for  the  eigenfrequencies  of  an  elastic  bar, 
i.e. 


cot  £  =  R£ 


(25) 


as  given  by  Timoshenko  (1955). 

For  R  =  0  (but  including  damping),  eq  20  and  21  reduce  to  those  given  by  Brown  and  Selway 
(1964),  i.e. 


f  -  \  Cte  -  1) 


and 


s  2sl"h'‘(ir) 

tan  -  = - 

2  n(2n  -  1) 


(26) 


(27) 


where  the  notation  has  been  changed  to  conform  to  this  report.  The  latter  two  equations  are  similar 
to  those  given  by  Lee  (1963)  if  the  small  angle  assumption  is  made  and  only  the  first  mode  is  con¬ 
sidered. 


COMPUTER  GENERATED  CURVES  AND  DISCUSSION 

Response  carves  Q  versus  ( for  three  modes  (Fig.  2-4) 

Figures  2-4  give  the  absolute  value  of  the  acceleration  ratio  of  the  free  end  of  the  bar  to  the 
driven  end  (defined  as  Q )  as  a  function  of  the  frequency  ratio  £  for  three  modes  of  vibration.  In 
generating  these  curves,  tan  S/2  was  arbitrarily  set  at  a  constant  value  for  each  curve  although  in 
a  real  material  one  would  expect  some  variation  with  frequency.  £  was  incremented  and  values  of 
Q  and  £  were  plotted  for  selected  values  of  mass  ratio  R.  In  each  mode  the  effect  of  end  mass  is 
obvious.  The  resonant  frequency  is  lowered  dramatically  with  increased  R  and  there  is  also  a 
decrease  in  Qmar  the  maximum  value  of  the  response. 

The  computer  plotter  was  programmed  to  print  a  plus  sign  on  the  Q  versus  curves  (see  Fig. 

2-4  and  Appendix  C)  when  the  phase  relationship  was  90° (corresponding  to  Re  =  0),  the  convenient 
experimental  point  discussed  in  the  section,  Measurement  ol  the  complex  modulus  at  90°  phase 
shitt  (p.  4).  This  point  corresponds  to  specific  values  of  and  Q'  also  discussed  previously.  Un¬ 
less  tan  5/2  is  very  small,  Q'  does  not  coincide  with  Qmax  but  is  shifted  to  the  higher  frequency 
side  of  resonance.  For  the  first  mode  with  R  =  0,  eq  20  reduces  to  cos  £'  =  0;  hence  £'  is  n/2  for 
any  value  of  tan  5/2  as  seen  in  Fig.  2.  However,  for  R  >  0  it  is  seen  that  £'  actually  in  '  -  eases 
with  increased  tan  5/2  for  a  given  mass  ratio  R  although  the  true  resonant  frequency  h,  lowered  with 
increased  tan  5/2.  Experimenters  must  not  confuse  Qmax  with  Q'. 

Variation  of  £'  with  R  for  various  values  of  tan  5/2  (Fig.  5) 

<f'  is  an  important  quantity  in  the  theory  since  it  is  used  to  compute  the  magnitude  of  the  complex 
modulus  E*  using  eq  16.  Figure  5  is  a  computer  generated  plot  of  the  variation  of  with  R  for 


Acceleration  Ratio 
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Figure  2.  Acceleration  ratio  Q  vs  frequen¬ 
cy  ratio  £  lor  various  mass  ratios  R  and 
values  of  damping,  tan  S/2,  first  mode. 


£,  Frequency  Ratio 


Figure  4.  Acceleration  ratio  Q  vs  frequency  ra¬ 
tio  £  for  various  mass  ratios  R  and  values  of 
damping,  tan  S/2,  third  mode. 


Figure  3.  Acceleration  ratio  Q  vs  frequen¬ 
cy  ratio  £  for  various  mass  ratios  R  and 
values  of  damping,  tan  S/2,  second  mode. 


Figure  5.  Variation  of  £  vs  mass  ratio  R  for 
various  values  of  damping,  tan  S/2. 
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various  values  of  tan  d/2.  decreases  significantly  with  increased  R  especially  for  the  higher 
modes.  For  any  given  value  of  R,  £'  increases  with  an  increase  of  tan  8/2  as  previously  seen  from 
Fig.  2.  Fig.  5  was  generated  from  eq  20;  tan  8/2  was  specified  and  for  an  incremented  value  of  R 
a  half  interval  search  method  (see  Kuo,  1965)  was  used  to  find  and  plot  the  values  of  £'.  The 
computer  program  is  given  in  Appendix  D. 

Tan  8/2  from  measured  values  of  Q'  for  typical  mass  ratios  (Fig.  6-8) 

These  curves  permit  the  determination  of  tan  8/2  from  experimentally  measured  values  of  <?' 
for  selected  values  of  mass  ratio  R.  The  computer  program  of  Appendix  A  used  to  solve  eq  20  and 
21  was  used  to  generate  these  curves. 

These  curves  plot  as  a  straight  line  on  log-log  paper  down  to  a  certain  value  of  Q'  (about  4.0 
for  the  first  three  modes)  and  then  show  a  slow  deviation.  This  implies  there  exists  a  relationship 
over  the  linear  range  in  each  mode  of  the  form 

Q'  tan  j?  =  i  (28) 

where  the  constant  C  is  a  different  number  for  each  value  of  R.  It  is  easy  to  show  from  Fig.  6  or 
directly  from  eq  20  and  21  that  for  A  =  0  the  constant  in  the  first  mode  is  n/2. 


Figure  6. 


Tan  8/2  vs  acceleration  ratio  Q\  lirst  mode. 
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EXPERIMENTAL  WORK 


The  experiment 

The  Objective  of  the  experimental  work  was  to  check  the  applicability  of  eq  20  and  21,  which 
include  the  mass  loading  effect  in  the  measurement  of  the  dynamic  modulus.  To  do  this  the  theory 
of  the  section,  Measurement  ol  the  complex  modulus  at  90°  phase  shift  (p.  4),  was  employed  using 
a  polymer  bar  to  simulate  a  soil  sample. 

The  experiment  consisted  of  driving  the  bar  with  a  small  end  mass  in  the  first  three  modes  and 
measuring  Q\  the  measured  acceleration  ratio,  and  /',  the  frequency  where  the  90°  phase  shift 
occurred.  The  test  was  then  repeated  with  a  larger  end  mass,  with  the  bar  shortened  to  maintain  a 
constant  frequency  in  each  mode.  The  computer  program  of  Appendix  A  was  then  used  to  compute 
E *  and  tan  S/2,  both  of  which  should  be  invariant  with  respect  to  mass  ratio  R  since  the  assumption 
of  the  theory  is  that  E*  and  tan  S/2  are  functions  only  of  frequency  and  are  independent  of  R.  The 
check  then  was  the  coincidence  of  E •  (or  tan  8/2 )  measured  at  various  mass  ratios  when  plotted  at 
the  constant  frequency. 

Apparatus  and  method 

A  schematic  diagram  of  the  testing  system  is  shown  in  Figure  9.  The  source  of  sinusoidal 
displacement  was  an  MB  Electronics  Model  EA  1500  electromagnetic  exciter  which  was  driven  by 
an  amplifier  and  audio  oscillator  (signal  generator).  A  frequency  counter  was  used  to  measure 
frequency.  Two  piezoelectric  accelerometers  were  employed  to  measure  acceleration  at  the  ends  of 


]  ACCELEROMETER 


TEST  BAR 


CHARBE 

AMPLIFIER 


VTVM 


VTVM 


PHASE  METER 


ACCELEROMETER 

\ - j 

MB  |__|  POWER 
EXCITER  [  (AMPLIFIER  I 


H  OSCILLOSCOPE  H 


FREQUENCY 

COUNTER 


SI0NAL 

SENERATOR 


Figure  9.  Schematic  of  the  testing  system. 
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Figure  10.  Test  specimen  assembly. 


the  bar.  The  output  from  the  accelerometers  was  fed  into  two  charge  amplifiers  and  read  from  two 
vacuum  tube  voltmeters.  The  amplified  signals  were  displayed  on  a  two-channel  oscilloscope  and 
the  phase  shift  was  observed  on  the  oscilloscope.  It  was  found  that  the  phase  of  the  two  signals 
could  be  accurately  measured  on  the  oscilloscope;  hence  in  the  latter  part  of  the  experiment  the 
phasemeter  employed  in  earlier  experiments  was  omitted. 

The  test  specimens  were  J/4-in.-diam  bars  of  low  density  polyethylene.  The  measured  specific 
gravity  of  this  material  was  0.915.  The  bars  were  used  as  received  from  the  supplier  with  no  heat 
treatment  and  only  the  ends  machined.  All  testing  was  done  at  75°F  ±  2°F.  The  method  of  fasten¬ 
ing  the  end  mass  and  fixing  the  bars  to  the  exciter  is  shown  in  Figure  10.  A  2-gram  accelerometer 
was  glued  with  Eastman  910  cement  directly  to  the  mass  or  in  the  case  of  the  low  mass  tests  directly 
to  the  free  end  of  the  bar.  Bar  lengths  are  given  in  Table  II. 

Calibration  of  the  accelerometers  was  done  at  the  three  frequency  ranges  of  interest  (data  were 
recorded  in  the  first  three  modes)  at  the  driving  magnitude  by  driving  the  accelerometers  back  to 
back.  Three  different  calibration  constants  were  used.  All  tests  were  made  with  the  exciter 
acceleration  set  at  10  G.  Phase  shift  in  the  electronics  was  carefully  checked.  The  transverse 
vibration  of  the  bar  was  checked  with  a  stroboscope  and  found  to  be  insignificant.  The  signals 
showed  no  visible  distortion. 

All  data  were  recorded  when  the  phase  angle  between  the  two  signals  was  90°.  The  value  of 
Q'  was  measured  and  the  frequency  was  recorded.  tan  S/2,  E* ,  Ej,  Eg  and  base  stress  were 
then  computed  using  the  computer  program  given  in  Appendix  A.  The  experimental  results  are 
presented  in  Figures  11  and  12  and  the  computed  results  are  presented  in  Table  II. 
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a.  End  mass  accounted  lor. 
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b.  End  mass  neglected  in  theory. 

Figure  11.  Experimental  results,  E*  vs  frequency. 


Frsquoncy,  Hz 

a.  End  mass  accounted  for. 


500  1000  5000 

Froqutney,  Hz 

b.  End  mass  neglected  in  theory. 

Figure  12.  Experimental  results,  tan  8/2  vs  frequency. 
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Table  n.  Computed  results  using  experimental  data.* 


Test 

no. 

Mass 

ratio 

R 

Bar 

length 

On.) 

Frequency 

r 

(Hz) 

Accel 

ratio 

Q' 

Frequency 

ratio 

e 

Tan  S/2 

E* 

fpsii 

E. 

(psi) 

E, 

(psi) 

i 

0.029 

10.02 

704 

8.21 

1.53 

0.0771 

72160 

71302 

11097 

2 

.029 

10.02 

2236 

2.53 

4.58 

.0811 

80843 

79780 

3 

.029 

10.02 

3842 

1.66 

7.64 

.0711 

85858 

84989 

12182 

35 

.029 

10.00 

712 

10.25 

1.53 

.0619 

73517 

72954 

9076 

36 

.029 

10.00 

2228 

3.40 

4.58 

.0610 

79952 

79357 

9737 

37 

.029 

10.00 

3830 

2.12 

7.64 

.0568 

84994 

84446 

9632 

38 

.029 

10.00 

3830 

2.18 

7.64 

.0553 

84995 

84475 

9384 

39 

.029 

10.00 

2215 

3.34 

4.58 

.0621 

79022 

78413 

9791 

40 

.029 

10.00 

707 

0.68 

1.53 

.0655 

72487 

71866 

9471 

41 

.029 

10.00 

718 

9.71 

1.53 

.0653 

74760 

74123 

9739 

42 

.029 

10.00 

2226 

3.09 

4.58 

.0670 

79807 

79092 

10659 

43 

.029 

10.00 

3839 

2.05 

7.64 

.0586 

85393 

84807 

9985 

44 

.029 

10.00 

702 

9.83 

1.53 

.0645 

71466 

70871 

9197 

45 

.029 

10.00 

2249 

3.06 

4.58 

.0676 

81460 

80716 

■  \m 

46 

.029 

10.00 

3882 

2.03 

7.64 

.0591 

87316 

86706 

■ 'or! 

59 

.489 

7.00 

722 

12.42 

1.08 

.0475 

73470 

73139 

6963 

61 

.414 

8.25 

2215 

2.43 

3.72 

.0530 

81364 

80908 

8600 

62 

.414 

8.25 

4067 

0.91 

6.85 

.0526 

86103 

85626 

9047 

79 

.504 

6.69 

738 

12.50 

1.08 

.0471 

71259 

70943 

6696 

85 

.410 

8.25 

3991 

0.99 

6.65 

.0489 

82887 

82491 

8089 

86 

.436 

7.76 

2326 

2.33 

3.70 

.0539 

80051 

79585 

8616 

*£*,  tan  5/2.  E*. 

E„  and  E2 

were  computed  using  the  computer  program  of  Appendix  A. 

Experimental  results 

The  ranges  of  E*  and  tan  S/2  which  were  computed  from  the  experimental  values  of  Q'  and  V 
for  various  mass  ratios  are  shown  in  the  upper  portion  of  Figures  11  and  12,  respectively,  at  the 
frequencies  corresponding  to  the  first  three  modes  of  vibration.  To  illustrate  the  errors  introduced 
in  computing  the  complex  modulus  ignoring  end-mass  effect,  the  lower  halves  of  these  figures  are 
plots  of  the  same  data  calculated  from  theory  neglecting  end-mass  effect.  All  data  are  taken  from 
Table  II  and  there  are  seven  points  at  each  frequency  (some  are  superimposed). 

It  may  be  concluded  that  large  errors  are  introduced  in  the  computation  of  E*  if  mass  loading 
effects  are  neglected.  For  example,  for  R  as  low  as  0.029,  E *  will  be  about  6%  low  in  the  first 
three  modes;  for  R  =  0.5,  E*  will  be  about  50%  of  its  true  value  in  the  first  mode.  This  is  illustra¬ 
ted  graphically  in  Figure  11. 

The  effect  of  mass  loading  on  the  computation  of  tan  S/2  is  seen  in  Figure  12.  The  errors 
introduced  using  R  =  0  theory  are  about  1%  for  R  =  0.029  but  become  greater  with  increased  mass 
ratio  and  mode  number.  For  example,  Young  (1967)  gives  an  error  of  130.6%  high  for  tan  S/2  for 
R  -  0.414  in  the  third  mode.  The  spread  in  the  experimental  data  for  tan  8/2  for  constant  values 
of  R  makes  it  difficult  to  interpret  the  data  in  these  tests. 


CONCLUSIONS  AND  SUMMARY 

The  theory  given  here,  including  end-mass  effect,  leads  to  more  nearly  correct  results  in 
computing  the  complex  modulus  from  vibrating  bar  test  data  and  should  be  adopted.  For  laboratories 
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using  the  90°  phase  shift  measurement  technique,  this  report  presents  curves  that  allow  direct  use 
of  experimental  data  to  calculate  the  complex  modulus.  For  experimental  data  that  fall  outside  the 
range  of  these  curves,  a  computer  program  is  presented  for  the  same  purpose. 
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iomc  (  i N/sr  c  ) . ?>x , 6nr  ( p:  i  i  .bx .  ;nf  1  .  i  x  , <;  p  ,£>•.; .gum  (x  r. , *.x , 

10HL  OF  MAR/) 

??  FORMAT  ( f- 1  0. 0  ,  ra  .?  .  2X  .  I  0  w’X  ,f  I  :  .  < ,  I  X  .  I  1  .  4  .  X  «  !  Vi'iXi  IT.'u,  |  /  ,  4  X  .  t  ’  . 

t  4X  ,f  1  0.6,'  X.TG.  ?  ) 

FN3 


c  typical  input 


C  R 

0 

OG 

L 

r 

N 

n  i  FPL 

0.0 

n.07 

0.011 

1  3.o 

0??C.o 

0  1 

i  .ov>r  -06 

0. 071  rv( 

n.o7 

0.P1  1 

C 

• 

r. 

o??r>  o 

0  1 

1  .oo-.f  -,v. 

o.o 

?0.6 

0.01  | 

13.o 

06030 

0? 

o . 07‘-.r  -"7 

0.07)03 

?r'.0 

0.P1  1 

.3.0 

06630 

0,0 

0  .  fl7or-0  7 

0.0 

17.0 

O.oi  1 

13.0 

1  0°M  3 

O'1 

0  .'.T'.f  -07 

0.071 03 

17.0 

o.°i  1 

1  3.  o 

l  oono 

07 

0  .  r.  3r  .f  -  0  7 

e  typical  Output 

R 

c 

FICPS) 

Vil/C 

TCEL/2 

C(  IN/SEC) 

0.0 

e.97 

2269 

1.571 

0.07C8 

126156 

0.072 

0.97 

2269 

1.466 

0.0705 

135140 

0.0 

25.60 

6639 

4.712 

0 .  COO  3 

123042 

0.072 

25.60 

6639 

4.4C7 

0.C079 

131567 

0.0  . 

17.60 

1C983 

7.054 

0.C072 

122170 

0.072 

17.60 

1C983 

7.369 

0.C065 

1  30160 

E  (  PS  1  I 

El 

E  2 

STRESS 

L  OF  OAR 

l 349692 

1336218 

190230 

2.691 124 

1  3.90 

1540047 

1 533538 

217226 

2.898381 

1  3.90 

1290243 

129C065 

21383 

0.980596 

1  3.90 

1475233 

1475047 

2  34  1  3 

1.099324 

1  3.90 

1271209 

1271076 

1  8381 

0.677436 

13.90 

1444042 

1443921 

18663 

0.815897 

1  3.90 
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APPENDIX  B.  STRESS,  STRAIN  AND  DISPLACEMENT  A8  A  FUNCTION 
OF  X  AND  R  (TABLE  I) 


COMPLEX  P,PL  .C2.GAMA,  XI  .  EPS  .  Gl  GMA  .  D  I  SPL  ,  H  1  C  .  OC  ,  yC  ,  EPS  I 
100  RE  AO (1*2)  Y.UZ.P,  AL.E  ,x .KOUNT 

?  format ( 3F io.o.f io.g.f  t  o .o ,f 1 o.f, 1 2 ) 

IE  (  KOUNT  )F>,  7.5 

5  WRITE (3,1 ) 

Y  =  2  •  *  A  T  AN ( Y ) 

AUX 

f)  1  =  x  *  T  A  N  (  Y/2.  ) 

R 1 Cs (0.0. I .0 ) *D1 
PL-A1-B1C 
P  =  PL-/AL 
GAMA  rP*PL 

C  Z-  ( CS 1 N<  PL ) +  G  A  M  A  *  C  C  0  5 (  PL  )  >  /  ( CCOS(PL ) -GA M A * C S I N ( PL  >  ) 
ZN=0.0 
11  A=ZN*X 

0  =  A*TAN( Y/2.  ) 

nc* in.o.i,n)«n 

XI  =  A  -UC 

fPSrU7»P«  (C2*CC0S(  XI  )  -C  S  I  N  (  XI  >> 
yc  =  <  o.o . i .0) *y 
f  PS  1  =F  PS*CEXP(  YC  ) 

SIGMArf;*LEXP<  YC)»EPS 

D I SPL=U7* (CCOG (  XI  !+C2*CGIN(  XI  >  -  1  .  > 

O I SPLS*CABS(D I SPL ) 

D  I  SPLR=REAl  (0  1  GPL ) 

D  I  GPL  I = A  I  MAG ( P  SPL) 

IE(DIC'PLR)  1G,1G,IS 
16  nnnAV=o.n 
GO  TO  2B 

1  6  DPHAf.r  =  AT  AN(  0  I  GPL  I  /O  I  SPLR  ) 

20  S I GMAS - CABS ( S I GMA  ) 

G I GMAP  sRC AL ( S I GMA ) 

G  I  G  M  A  1  r  A  I  V A G ( S I GMA ) 

1  E  (  S I GMAP  )  17, IB, 17 

IV  OPHAf-r  -  AT  AN  (  G  I  GMA  I  /SI  OMAR  ) 

IB  EPSG  =  C  APS ( EPG ) 

EP  >R  =  PE  AL ( EPS ) 

EPS  I  = A  1  MAG (EPS ) 

IE(EPGR)  10,20,10 
1  o  rpMftjr  r  AT  AN( !  Pc  I  /f  PGP  I 
20  W  P I T  T (3,3)  PI  GPL  ,G  I  gma ,>  pg j 

y.’R  IT  f  (  3 ,2  1  )  r  I  GPL  ,  g  i  r.M  ag  ,  <-  l’G^ 

WR  1  TE  (  3 ,2  1  )  PPM  A  '■  f  ,Gr,,iAGr- ,ppn,'GP 

7N-7N+0.7‘- 

IE ( 7N- 1.20)  11,17,12 

1  2  e;R  Iff  (  3,(i ) 

GO  TO  100 
7  C All  EXIT 

1  FORMAT  (  l  H  1  ,GX  .  1  2H0  I  GPL  AC  *  ME  NT  ,  poy  ,6  IP  T  R-  GG  ,  2  3X  ,  6M  G  T  R  A  l  N  /  ) 
3  FORMAT  (  /2f,  1  2  .  •• ,  6X  .  2G 12.'.,  E  X  .  2 ' 1  2  .  A  ) 

2  1  (-  OPM  AT  (  o  1  2  .  A  ,  1  2X  ,  0  1  2.  A  .  >  2X  ,  o  1  7  ,  A  t 

6  FORMAT(//) 

A  F  PPM A T  (  I G  ) 

END 


OUTPUT  FPS  THIS  DAI 
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APPENDIX  B 
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APPENDIX  C.  Q  VERSUS  £  FOR  VARIOUS  R  AND  TAN  8/ 2  VALUES 

(FIG.  24) 


< 

a 

G 

o 

a 

L 


O 

OOOOOCO  OOG'CCOO 
X  0~(\,cc~<\.  O-CJCO-M 

— •  •••••••  ••••••« 

c.  ocooooo  ooooooo 


u 

a 

a 

< 


a 

O  OOOOOOO  CD  CO  <D  CD  CD  Oj  CD 

u.  •••••••  ••••••• 

nnnrinnn  tn  40  10  tr>  in  «f)  in 

< 


< 

c 


00©  —  inirinr*ooo  —  intnirr* 

•  ••••••••••••••• 

000000000000000000 


* 


CALL  PLOT  C  1 00 ) 

CALL  FX1T 

FORMAT (  1H1  *5X.  1  HP,  6X»3HTANDEL/2*,3X«  1HX«9X,  l  HO/  ) 
FORMAT ( AF 1 0.5l 


APPENDIX  D.  VERSUS  R  FOR  VARIOUS  VALUES  OF  TAN  8/2  (FIG.  5) 


1  C  =  1 

C 

DATA  fOP 

APPENDIX  1) 

PROGRAM 

tCC=90 

C 

DO  11  1=1,2 

Y 

P 

X 

READ! 1,22) 

1 1 

WRITE(3,22) 

0.8 

l  •  n 

0.0 

01 

22 

FORMAT! 1X.79H 

0.8 

i  .n 

0.  1 

0  1 

1  > 

0.8 

i.n 

0.2 

01 

EPS=0.00l 

3.4 

b.0 

0.0 

or* 

CALL  PLOT! |C. 0.0, 0.6. 6. 0,0, 1 .0.0,10.0,10,0.1.0) 

3.4 

*>•0 

0.  1 

or 

CALL  PLOT (99) 

3.4 

r».o 

0.2 

or 

30 

READ! 1.3)  XI 1 .X22.Y.N 

6.5 

9.0 

0.0 

o:< 

IF  I  X 1  1  )  6.7,6 

6.6 

9.0 

0. 1 

0T 

6 

WR1TEI3.4)  N 

6.6 

9.0 

0.2 

0T 

R=-0. 00125 

O 

• 

o 

DO  25  K  =  1 .6 
DO  10  1*1,90 
X  1  *  X  I  1 
X2=X22 
R=R+O.OOI 25 

PI =3. 14 159 

COT  1 =  C06 ( X  I  )/5IN(XI  1 
C0T2=C05(X?1/SIN(X?  ) 

FX  1  =  COT 1-R»X1 +  R*Xl «Y»C0T1 »TANH(Xl *V ) 
FX2=C0T2-R*X2+R«X2«Y»C0T2*TANH(X2*Y ) 
70  X=(Xl+X?)/2. 

CO  X=COS(XI/S!N(x > 

FX--C0TX-R*X  +  R«X*  Y*COTX«TANH(X*Y  1 
IF(ADS(FX)-EPr. )  15,15,60 

60  IF(FX*FX1>  50,16,80 
60  X2  =  X 
Fx2=rx 
GO  TO  70 
BO  XI =X 

fx i -rx 

GO  TO  70 

15  CALL  PLOT! ICC, R, XI 
10  CONTINUE 

WRITE (3,1 )  Y , R , X 
25  CONTINUE 

CALL  PLOT (99) 

GO  TO  30 
7  CALL  PL. OT  (  100) 

CALL  EXIT 
1  FORMAT (5r 10.5) 

3  FORMAT! 3F10.5, 12) 

4  FORMAT ( //6X.4HN  =  ,1?//) 

END 
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